Sharp and Infinite Boundaries in the Path Integral Formalism 



We revisit the analysis of sharp infinite potentials within the path integral formalism using the image 
method pQ. We show that the use of a complete set of energy eigenstates that satisfy the boundary 
conditions of an infinite wall precisely generates the propagator proposed in Ref. pQ. We then show the 
validity of the image method by using supersymmetric quantum mechanics to relate a potential without 
a sharp boundary to the infinite square well and derive its propagator with an infinite number of image 
charges. Finally, we show that the image method readily generates the propagator for the half-harmonic 
oscillator, a potential that has a sharp infinite boundary at the origin and a quadratic potential in the 
allowed region, and leads to the well known eigenvalues and eigenfunctions. 

PACS numbers: 03.65.-w 

Introduction 

Sharp and infinite potentials are difficult to handle in the path integral formalism [21 El E] - Goodman [T] , 
in one of the first papers to address this problem, advocated the use of the image-point method to account 
for paths that wander into the forbidden region. He applied this intuitive method to the infinite square well 
and generated the correct set of eigenvalues and eigenfunctions. In Ref. [5], Sokmen provided additional 
support for the image method by connecting the infinite square well (ISW) to the Rosen-Morse potential via 
a change of variables (point canonical transformation). Starting with the Rosen-Morse potential, which has 
no sharp boundary, Sokmen derived the propagator of Ref. [T], and thus verified, at least for the case of 
infinite square well, the validity of the image method. Auerbach and Kivelson [B], Nevels, Wu and Huang [7], 
and Auerbach and Schulman [8] have provided much more rigorous reasoning to support the image method 
and have extended it to systems with non-zero potentials in the allowed region . 

In this paper, we aim to provide further support for the validity of Goodman's conjecture. First, we show 
that the use of eigenfunctions that vanish at the infinite boundary automatically generates the propagator 
Goodman found for a particle in a restricted domain, i.e. it naturally generates the propagator with an 
image point. We then, as Sokmen did in Ref. [5], generate the correct propagator for the ISW that contains 
an infinite number of image charges, starting from the Rosen-Morse potential. However, instead of using a 
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Figure 1: Inclusion of the term originating from the image-point (— x, t) has the effect of excluding all 
forbidden paths such as the one shown by the solid line. 

point canonical transformation [9l [10l [11] to map the Rosen-Morse potential to the ISW, we show that a 
simple change of parameter is sufficient for this mapping. We use supersymmetric quantum mechanics to 
show the validity of this limiting procedure. Finally, we use the image method to derive the propagator for 
a half-harmonic oscillator with a sharp boundary at the origin and show that it generates the well known 
spectrum for this system. This final case provides an explicit example of a particle in a non-zero potential 
with a sharp boundary. 



The propagator for a free particle in a restricted domain 

In this section we will derive the propagator for a particle in the presence of an infinite potential barrier using 
the path integral approach. The particle is free in the allowed region: < x < oo. The Hamiltonian for this 

~ - 2 

system can be written as H = |^ + V(x), where V(x) is given by 

!oo for x < 
" ■ (1) 
for x > 

Our goal is to determine the propagator K^{x' 1 1'\ x, t) = u(x' ,t'\ x, t)-% = {x'\ e~^ H ^ ~^\x) for this system 
that respects its boundary conditions. Wc have used the notation | x, t)u to denote an eigenstate of the 
position operator in the Heisenberg representation that is related to an eigenstate in the Schrodinger picture 
by | x, t)fi = e lHt / h \ x). To simplify the notation, we will drop the subscript %. 
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Following Rcf. we split the time-interval (t'—t) into n+1 segments, each of size e, letting denote the time 
t + £e, with t + (n + 1) e = t' . Now, using the completeness of state vectors | ag, fc), i.e., J da^ | (ajj, fe| = 1, 
where denotes position at an intermediate time we can write the matrix element as 

(x',t'\x,t) = Jdx n [dx n -! J- ■ ■ J dx x (x',t'\ x n ,t n ) ■ ■ ■ (xt,tt\ xt-\,U-\) ■ ■ ■ {xx,t\\ x,t) . (2) 

For very large n, and hence for an infinitesimally small e, the propagator (xg, t(\ a^-i, tg_i) can be written as 

dE (x e \e-^\E) {E\x t -i) 



(xe,te\ xg-utt-i) = (ass 







e » 









d£ (^l-E) (^|a;/_i) e- ie ^ 



2dfe -ie^&± 

sinto smA^a;^_i e h 



(3) 



where the function H(pi), without a hat, represents the classical Hamiltonian J^-. The kets represent 
the energy eigenstates for a free particle in the domain < x < oo. The corresponding coordinate repre- 
sentation is given by (a;| E) — y ^ k sinfcx, which vanishes at the infinite wall at x = 0. The wavenumber 
k is related to the energy by E — 2m . The factor in front of sin kx is derived from the requirement 
that J dE(x\ E)(E\ y) must equal 5(x — y), the completeness condition. To compute the above infinitesimal 
propagator {xe, te | xi-i, ti-i), following Ref. [7J, we write J °° sin^aj sin£^a;£_i as 



sin sm^in) 



dh 



g i*?(a(-x/_i) _ ife(:i*+x*_i) 



Substituting this expression into Eq. ([3|) yields 

oo 2-7T 



hi.pi) 



Inserting this infinitesimal propagator into Eq. (2), we get 

n ,.oo n+l r „oo 



^ /"OO ~i~ |>OG 11 

{ x>,t>\x,t)=ii d „n / S 

These Gaussian integrals can be carried out using the identity dx e 

n /.oo n-\-l 1/2 



. They yield 



^im(x£—xi-i) /2he ^im(x£-\-X£-i) /2he 



(4) 



(5) 



Terms in the square bracket, upon integration, lead to free-particle propagators from (x,t) (x',t') and 



3 Such a substitution is justified because the transformation k — > —k has no effect on the Hamiltonian as it is a quadratic 
function of the momentum. 
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(— x, t) — > {x' ,t') respectively I.e. 

/ m x 



(x',t'\x,t) = 



\2mh(t' - t) 



e im(x' -xy /2h(t' -t) _ e im(x' +xY /2h(t' -t) 



={x', t'\x,t) F - (x', t'\-x,t) F = K F (x', t'\ x, t) - K F (x', t'\ -x, t) , 

where the subscript F denotes a propagator for a particle that is free in the region — oo < x < oo. Thus, we 
find that the propagator Kr(x' , t'; x, t) in the restricted domain < x < oo is given by the difference of free 
particle propagators originating from the point {x,t) and its image point (— x,t), as suggested in Ref. [JJ. 



Green's function for the ISW from the Rosen-Morse potential 

Supersymmetric quantum mechanics (SUSYQM) |14j [15j [16j H3 HH] is a generalization of the ladder operator 
formalism for the harmonic oscillator. It makes use of first order differential operators = (=F ^ + W(x, a)) , 
where the superpotential W(x, a) is function of x and a parameter a, and embodies the interaction in this 
formalism (instead of the potential.) The operators A are hermitian conjugates of each other and help us 
define two partner hamiltonians 

h - = A±Ar = + H + "<*■•>) - + * ^ - + v ^> • 

where we have set h = 2m — 1. A superpotential W(x,a) generates two Hamiltonians H± with potentials 
V±(x, a) given by V±(x) = W 2 (x,a) ± ^f. For example, the superpotential W(x) = \ ujx generates the 
partner potentials V± (x) = | ujx 2 ± |u; both describe one-dimensional harmonic oscillatorfj^J 

The infinite square well is an extensively studied system in SUSYQM 19 . The general Rosen-Morse and 
Poschl- Teller potentials reduce to ISW for a particular value of a parameter. We will use this relation with 
the Rosen-Morse potential to derive the propagator for the ISW. We begin with the superpotential for the 
Rosen- Morse system: W(x, b) = —bcotx. The domain of this potential is < x < ir, and the partner 
potentials V±{x,b) are given by 

V±(x, b) = W 2 (x, b) ± 4- W(x, b) = b(b ± 1) cosec 2 x - b 2 . (6) 
da; 

For an arbitrary positive value of 6, the supersymmetry is said to be unbroken [TBI [T71 [TB] , and the potential 
V- (x, b) — b{b — 1) cosec 2 x — b 2 holds a zero energy ground state with corresponding eigenenergies given by 
En ^ = (b + n) 2 — b 2 . Note, that for 6 = 1, this potential reduces to V-(x, 1) = — 1. Thus, V-(x, 1) describes 
an Infinite Square Well of width n, and with its bottom at —1. We will show that for limfj-j-i the Green's 
function for V— (x, b) reduces to that of an ISW with an infinite number of image charges. 



4 See pages 31-37 of Ref. g]. 

5 Note that for harmonic oscillator, both partners, V+(x) and V_(x), have the same x-dependence. This feature, known as 
shape invariance, plays an important role in the exact solvability of this system [21 22 . 
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In Refs. [S1Q2], the authors derived the Green's function G(xf, X,*, E) for the potential Vps(x, s) = (s 2 — |j cosec 2 
It was given by 

G(x f ,x i ,E) = Y j ( n + s + 2 \ T ^ 2 ^ 1 ' ) (ana/)* P^cosx/) ((sinx,) 1 P^ (cos x,)) * , 

where represents the Associated Legendre Polynomials. For s = 6—5, the potential Vps(x,s) and 
V~(x,b) differ only by the constant — b 2 , and hence, they both must have identical Green's functions, the 
only difference being that the eigenvalues will differ by that same constant. For 6=1, i.e. s = |, the Green's 
function is given by 

G(x f ,Xi,E) = V „ r ( n + 2 j (sina?/)' P~K (cos a;*) ((sinx,)' P~* (cosx t ))* . (7) 

Now, we use the identitjj^] 

P^) = ^—X ( Z ^Y F (-u,t, + (8) 

to get 



r(i-/x) V^-! 



D -i , 1 /cOSX/ + l\ * / 1 3 3.20;/ 

P 2 ! (cosxf) = — t^v F \ -n , n + -; -; sin — 

«+2 v /; r(|) Icosxf-i/ I 2' 2' 2' 2 



■■.2J 



- 1 c ** W Sln ^ sin[(n+l)ay] 

A V cos ^ (n+l)sin(^) ' 1 j 

In deriving above, we have used the identity 7 F (a, 1 — a; |; sin 2 z) = ^^""^j ■ Substituting the above 

_ A _ A f~ 

expression for P n+ 2 i (cosxy), and a similar expression for P n+ 2 i (cosx^), into Eq. (7), we get 



00 (n + 1) 2 2 

G(x f ,Xi,E) = V -2 sin[(n + l)x/] sin [(rc + l)xi] 

„=o ^ - (n + 1) l n + tj 7T 



(10) 



J5 n i'nixf) lp n (xi) 

which clearly produces the correct spectrum for the infinite square well with ground state energy equal to 
6 2 = 1. 

Thus, we were able to connect the Rosen-Morse potential to the ISW simply by considering the system for 
a specific value of the parameter b. By taking the limit of b — > 1, we have avoided employing the point 
canonical transformation (PCT) of Ref. [5], which required changes of both the dependent and independent 
variables This much simpler method of mapping the Green's function of a continuous potential to the 

6 Abramowitz & Stegun - pp. 332, eq. (8.1.2). 
7 Abramowitz & Stegun - pp. 556, eq. (15. 1.15). 

8 It is worth pointing out that point canonical transformation often provides a way to derive information about one system 
by connecting it with another system 1201 . A substantive research on point canonical transformation in path integral formalism 
was carried out in Refs. l9l 1101 [TT1. 
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Green's function for a system with two sharp, infinite boundaries has allowed us to generate the propagator 
containing infinitely many image points with considerable ease. 



The propagator for a half-harmonic oscillator potential 

As a further verification of the image method, we will determine the propagator for a half-harmonic oscillator 
(H-HO) using the image method, and thence its spectrum. In Ref. [5], the authors use the path decomposition 
expansion (PDX) [6] to validate the image method for a general potential in the allowed region. We seek 
to verify their proposition by showing that the image method readily produces the propagator for the half 
harmonic oscillator, a potential defined over the half line < x < oo. Only the odd energy levels of the 
harmonic oscillator are expected to correspond to the eigenvalues and eigenfunctions of the H-HO, since the 
even levels would not satisfy the boundary condition that the eigenfunction must vanish at x = 0. 



The H-HO can be viewed as a system defined by 

V(x) = 



hkx 2 



for x < 
for x > 



containing a sharp and infinite boundary at x = 0. 



We start with the well known propagator Ky lo (xf,tf,Xi, ii^jjfor a harmonic oscillator defined over the entire 
real line (— oo < x < oo). It is 

imuj 



K ho (xf,tf,Xi,ti) 



2irhi sinwr 
where we have substituted r for (it — U). 



exp 



2h sine 



[(x^+x?) cos lot — 2xi xj\ I, 



(11) 



Then, assuming the applicability of the image method, the propagator for the H-HO is 
fCh-ho (xf,tf,Xi,ti) = Kho {xf,tf,Xi,ti) - K ho (x f ,tf, -Xi,U) 



2nhi sinwr V 2h sinwr 



ex P I 77T-J—Z [{ x f + 



2 1 X?) COSWTl I X 



exp 



2 imuj XiXf 



exp 



2 imu XiXf 



2h sincjr J \ 2H sinwr 

We will now show that this expression generates the correct eigenvalues for the H-HO. We start with 



K (x, tf,X, U) = 2J e R En T 4>n{xf)4> n {Xi) 



(12) 



(13) 



From here, we substitute This substitution and a subsequent integration allows us to identify 

the eigenvalues: 



/>00 00 /"OO 00 

/ dxJC(x,t f ;x,ti) = J2e~ iEnT / dx \M%)\ 2 = J2 e ~ iEnT ' 

Jo ,i=0 Jo n=0 



(14) 



3 See page 200 of Ref. 0. 
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where we have used the orthonormality of ip n (x); i.e., J dx |^f n (a;)| = 1- Thus, the propagator Kh-ho ( x itf) x i U, 
is then given by 

„2 



Kh-ho (x,t f ;x,ti) 



itilo 



2nhi sine 



h sine 



imuix \ I imoox 
exp ( 1 —. (coscjt — 1) — exp - — . (coswr + 1) 



h sine 



Then, as in Eq. (14 1, we integrate the Gaussian Kh-ho (x,tf,x,U) 

1 



da; Kh-ho (x,tf,x,ti) 



Further simplification leads to 



2 V 2-nhi sine 



7rz h sinwr 



7rz h sin LOT 



dx Kh-ho (x,tf,x,ti) 



moo (cos lot — 1) y mio (coswr — 1) 



1 1 

y/2 (COS WT — 1) y/2 (COS LOT + 1) 



(15) 



(16) 



Now, using 1 — cos wt = 2sin 2 ( < rf ) and 1 + coswt = 2cos 2 {^&\ we get 



da; Kh-ho (x,tf,x,ti) 



2i sin *f 2 cos 



1 2e" — 



2e 



_ 3im 

e 2 



2CJT c —%(jJT 



1-e" 



_ 3to; 

e 2 



n=0 



Thus, we see that we get the expected eigenvalues for the half-harmonic oscillator 



E n = [2n+- )ftw, n = 0, 1,2,3, • 



(17) 



Thus, the energy levels of the H-HO emerge to be the same as the odd energy levels of the full harmonic 
oscillator. 



It is worth pointing out that our result could be arrived at much easier if we had started with 
K ho (x f ,t f ;x u U) =jr / e-*( n+ i> T H n (J^fxf) e"^ x ^n 

n=0 ^ " ' 

The propagator for the H-HO is then given by 

Kh-ho (xf,tf,Xi,ti) = K ho (xf,tf,Xi,ti) - Kho (xf,tf, -Xi,ti) 



(18) 



n=0 



-i(n+|)wT -tffi (x 2 f +xf) 



H n [ \l -z- Xf ) H n [ \ / — X,, 



//, ( \j — x f )H n 



(19) 
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Since H n {—x) = (—l) n H n (x), we find 



Kh-hoix^tf ;x t ,U) = J2^ in+i)ulT ^ W{x2f+xl) [l-(-l) n ]fT» (vT^J H » V V /, 

n odd ^ V / \ V / 



n=0 



muj 



<t>n(xf) <t>n(xi) 

(20) 

Note the factor of y2 in the normalization is necessitated by the change in domain. Thus, by comparing this 



result to Eq. ( 13 ), we see that the energy eigenvalues and corresponding eigenfunctions for the half-harmonic 



oscillator are given by 



E n = ^2n+ ^ Hu and 4> n (x) = V2e ™ x * H 2n +i (J^J^ x ^ 



respectively. 



Thus, beginning with two different initial points, we have shown that the image method readily generates 
the correct spectrum for the half-harmonic oscillator. 



Conclusion: 

We derived the propagator of a free particle in a restricted domain (0 < x < oo) using eigenstates that 
vanish at the infinite boundary, showing that image contributions naturally emerge from this method. We 
also showed that the Rosen-Morse potential, for a specific value of a parameter, reduces to the Infinite 
Square Well potential, and that the corresponding Green's functions map into each other. Thus, starting 
from a well defined potential with no sharp boundary, we derived the Green's function for a system with 
two sharp boundaries and an infinite number of image charges. The power of our method lies in the fact 
that by employing a change of parameters there was no need for a point canonical transformation as in Ref. 
[5 . Finally, we used the image method to derive the propagator for a half-harmonic oscillator and showed 
its correctness by generating a correct set of well known spectra for this system. This final case provides 
an explicit example of a particle in a non-zero potential with a sharp boundary and supports Aurbach and 
Schulman's proposition in Ref. [8 . 
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